Cooper Pairs' Magnetic Moment in MCFL Color Superconductivity 
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We investigate the effect of the alignment of the magnetic moments of Cooper pairs of charged 
quarks that form at high density in three- flavor quark matter. The high density phase of this matter 
in the presence of a magnetic field is known to be the Magnetic Color-Flavor-Locked (MCFL) phase 
of color superconductivity. We derive the Fierz identities of the theory and show how the explicit 
breaking of the rotational symmetry by the uniform magnetic field opens new channels of interactions 
and allows the formation of a new diquark condensate. The new order parameter is a spin-1 diquark 
condensate proportional to the component in the field direction of the average magnetic moment of 
the pairs of charged quarks. In the region of large fields, the new condensate's magnitude becomes 
comparable to the larger of the two scalar gaps. Since there is no solution of the gap equations with 
nonzero scalar gaps and zero value of this magnetic moment condensate, its presence in the MCFL 
phase is unavoidable. This is consistent with the fact that the extra condensate does not break 
any symmetry that has not already been broken by the known MCFL gaps. The spin-1 condensate 
enhances the condensation energy of pairs formed by charged quarks and the magnetization of 
the system. We discuss the possible consequences of the new order parameter on the issue of the 
chromomagnetic instability that appears in color superconductivity at moderate density. 

PACS numbers: 12.38.Aw, 12.38.-t, 24.85. +p 
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High dense quark matter 



o 

=5 



(N 
> 
00 
&\ 

(N 

in 
o 



I. INTRODUCTION 

It has been suggested long ago that a phase transition 
to color superconducting (CS) quark matter may take 
place at high baryon density and sufficiently low temper- 
ature [l|. In nature, the most likely place to accomplish 
CS is the core of compact stars. The typical density in- 
side a compact star can be of the order of several times 
the nuclear density. In such media, quarks may be re- 
leased from hadrons and thanks to the low star's temper- 
ature caused by the rapid cooling through the neutrino 
emission taking place after the supernova collapse, they 
form CS Cooper pairs. 

Compact stars, on the other hand, are very magnetized 
objects. From the measured periods and spin down of 
soft-gamma repeaters (SGR) and anomalous X-ray pul- 
sars (AXP), as well as the observed X-ray luminosities 
of AXP, it has been found that a certain class of neu- 
tron stars named magnetars can have surface magnetic 
fields as large as 10 14 — 10 16 G |2|. Moreover, since the 
stellar medium has a very high electric conductivity, the 
magnetic flux should be conserved. Hence, it is expected 
an increase of the magnetic field strength with increasing 
matter density, and consequently a much stronger mag- 
netic field in the stars' core. Nevertheless, the interior 
magnetic fields of neutron stars are not directly accessible 
to observation, so one can only estimate their values with 
heuristic methods. Estimates based on macroscopic and 
microscopic analysis, for nuclear Q, and quark matter 
considering both gravitationally bound and self-bound 
stars [H, have led to maximum fields within the range 
10 18 — 10 20 G, depending if the inner medium is formed 
by neutrons or quarks 



The presence of a strong field in the star core can in 
principle modify the properties of the matter phase there 
and lead to observable signatures. Therefore, the inves- 
tigation of the properties of very dense matter in the 
presence of strong magnetic fields is of interest not just 
from a fundamental point of view, but it is also closely 
connected to the physics of strongly magnetized neutron 
stars. 

The effect of a strong magnetic field in the phase struc- 
ture and gap magnitude of a color superconductor was 
first studied in Q. There, it was shown that the pres- 
ence of a magnetic field changes the color-flavor-locking 
(CFL) phase that realizes at large densities, producing a 
difference between the gap that gets contributions from 
pairs of oppositely charged quarks, denoted in Q by As, 
and the one that only gets contributions from pairs of 
neutral quarks, denoted by A. Charged and neutral in 
this context refers to the electric charge associated to the 
"rotated" or in-medium electromagnetic group that re- 
mains unbroken in the CFL phase [6|. In this phase, even 
though the original electromagnetic U(l) em symmetry is 
broken by the pairing, there is a residual symmetry that 
gives rise to a new long-range field composed by a linear 
combination of the conventional photon field and the 8th 
gluon field [y, [7| ■ The CFL pairs remain neutral with re- 
spect to this new "rotated" electromagnetism, and hence, 
the color superconductor can be penetrated by a "ro- 
tated" magnetic field without the restriction of the Meiss- 
ner effect [8j . The mixing angle between the original elec- 
tromagnetic field and the 8th gluon is such that a regular 
magnetic field penetrates the superconductor through its 
rotated magnetic field component with almost unchanged 
strength. The new phase that forms in the presence of 
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the magnetic field also has color-flavor-locking, but with 
a smaller symmetry group SU(2)c+l+r an d this change 
is reflected in the splitting of the gaps. This magnetic 
phase was called Magnetic CFL (MCFL) [|. Similar to 
the CFL case, the MCFL phase is invariant under the 
rotated electromagnetic group. 

At sufficiently strong magnetic fields (eB > /j, 2 ), the 
separation between the two gaps of the MCFL phase be- 
comes significant, as the energy gap with contributions 
from pairs formed by rotated-charged quarks raises with 
the field while the other gap decreases with it [5(. How- 
ever, contrary to what our naive intuition might indicate, 
a magnetic field does not need to be of the order of the 
baryon chemical potential to produce a noticeable effect 
in a color superconductor. As shown in the color su- 
perconductor is characterized by various scales and dif- 
ferent physics can occur at field strengths comparable 
to each of these scales. Specifically, the MCFL phase is 
entirely attained when the field strength induces masses 
equal or larger than those of the charged Goldstones' con- 
stituents (i.e. twice of the CFL gap), and hence allowing 
the decay of this bound states. This effect is manifested 
at moderate magnetic fields (henceforth magnetic field 
always means rotated magnetic field), through the fact 
that all the energy gaps exhibit oscillations with respect 
to eB/fi 2 < 1 jlfllllj]. This behavior is due to the well 
known de Haas-van Alphen effect produced by the dis- 
creteness of the Landau levels in an external magnetic 
field, which becomes noticeable at field values beyond the 
threshold field needed to induce the Goldstone-mesons 
decay. 

Although there are several interesting alternatives for 
CS pairing patterns possessing nontrivial magnetic prop- 
erties as well [l2j, in this paper we restrict ourselves to 
the color antitriplet channel that is known to be favored 
for three massless quark flavors at high densities. In the 
presence of a magnetic field the ground state of this sys- 
tem is in the MCFL phase 0. However, it should be 
noticed that in this phase the Cooper pairs formed by 
charged quarks have nonzero magnetic moment because 
the quarks in the pair not only have opposite charge but 
also opposite spin. It is reasonable to expect that the 
MCFL ground state could have a net nonzero magnetic 
moment Jl which could couple with the applied magnetic 
field as ~ Jl- B. This extra magnetic energy contribution 
of the diquark ground state would be reflected in the ex- 
istence of a new order parameter Am in addition to A 
and A b • The main goal of the present paper is to explore 
precisely this possibility. 

Symmetry arguments can give additional insight on the 
likelihood of an extra gap Am in the MCFL phase. The 
presence of a magnetic field breaks the spatial rotational 
symmetry 0(3) to the subgroup of rotations 0{2) about 
the axis parallel to the field. As will be shown in this 
work, this fact opens new attractive pairing channels that 
are not available in the CFL phase. One of these channels 
has symmetric Dirac structure C757172. As it will be dis- 
cussed in detail later, a condensate with this structure is 



a spin-1 condensate with zero spin projection in the di- 
rection of the magnetic field. A diquark condensate with 
this new structure would have the same energy dimension 
as the regular gaps and would account for a magnetic en- 
ergy associated to the projection of the magnetic moment 
of the ground state along the external field. With this 
understanding, we will call it the magnetic-moment con- 
densate, because it should be proportional to the mag- 
netic moment of the ground state. We will see that this 
interpretation is in agreement with the way Am enters in 
the dispersion relations of the charged quasiparticles. As 
it will be shown, the color-flavor structure of Am can be 
chosen in such a way that this condensate does not break 
any symmetry that has not already been broken by the 
gaps A and A^. Hence, such a magnetic-moment con- 
densate is not symmetry-protected and in principle can 
be different from zero. This resembles the reasons that 
led to the existence of the symmetric gaps in the CFL 
phase. The difference is that the color channel for the 
symmetric gaps is repulsive, so their magnitudes are typ- 
ically much smaller than that of the antisymmetric ones. 
This is not the case with the magnetic moment, because 
the color channel here is going to be the same antitriplet 
attractive channel that favors the antisymmetric gaps. 
Nevertheless, since the magnetic-moment condensate is 
a direct consequence of the external magnetic field, we 
will see that its magnitude becomes comparable to the 
energy gaps only at strong field values. 

Despite some fundamental differences, the MCFL sce- 
nario described above shares a few similarities with the 
dynamical generation of an anomalous magnetic moment 
recently found in massless QED [Kjl . Akin to the Cooper 
pairs of oppositely charged quarks in the MCFL phase, 
the fcrmion and antifermion that pair in massless QED 
also have opposite charges and spins and hence carries a 
net magnetic moment. A dynamical magnetic moment 
term in the QED Lagrangian does not break any symme- 
try that has not already been broken by the chiral con- 
densate. Therefore, once the chiral condensate is formed 
due to the magnetic catalysis of chiral symmetry breaking 
EMU, the simultaneous formation of a dynamical mass 
and a dynamical magnetic moment is unavoidable [l3| . 
As we shall show below, the simultaneous generation of 
the MCFL gaps and the magnetic-moment condensate 
Am is also unavoidable in the MCFL phase of color su- 
perconductivity. There is simply no consistent solution 
with nonzero energy gaps and zero Am- 

More important, even though the solution for Am is 
negligible at small magnetic fields, it becomes compara- 
ble to Ab at strong fields, where the LLL contribution 
is more relevant, indicating that the main contribution 
to Am comes from the pairs formed by quarks that lie 
in the LLL. In addition, the gap As is actually larger 
(about twice at very strong fields) when Am is taken 
into consideration in the analysis than when it is ignored. 
These two facts make the new parameter Am of partic- 
ular relevance for the potential realization of the MCFL 
phase in the core of magnetars. As known, the CFL 
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phase may become chromomagnetic unstable at interme- 
diate densities due to the pairing stress generated by the 
strange quark mass and the color and electric neutral- 
ity conditions jl7j |. In this context, any effect that can 
augment the effective gap magnitude will contribute to 
stabilize the phase by pushing the emergence of the in- 
stability to smaller regions of densities. Since magnetars 
have the strongest surface fields, they should also have 
the strongest fields in the core, so they are the best can- 
didates for the realization of the magnetic CFL phase 
studied in this paper. The larger value of the gap at 
strong fields will be also reflected in a larger critical tem- 
perature for this color superconducting phase, a property 
that opens yet another possibility for the realization of 
the MCFL state because the conditions of high densities, 
low temperatures, and strong magnetic fields will likely 
coexist in the planned low-temperature/high-density ex- 
periments at NICA@JNIR, CBM@FAIR and low-energy 
RIHC [3. 

In this paper, we investigate the generation of a mag- 
netic moment condensate in the MCFL phase of color 
superconductivity, that is, in a three-flavor quark sys- 
tem at high density and in the presence of an external 
magnetic field. We perform our calculations within a 
Nambu-Jona-Lasinio (NJL) model in the presence of an 
external magnetic field. The outline of the paper is as 
follows. In Sec. II, we discuss the new pairing chan- 
nel that becomes available in the MCFL phase due to 
the external magnetic field, allowing for a new param- 
eter (the magnetic-moment condensate Am) in the gap 
matrix. This new gap is a spin-1 condensate with zero 
spin projection (Ms = 0) along the field. To explore the 
simultaneous solution for A a/, A and As, we introduce 
a new MCFL ansatz for the color-flavor-Dirac structure 
of the gap matrix that contains the two usual energy 
gaps and the magnetic-moment condensate. This matrix 
respects all the symmetries of the old MCFL phase Q, 
i.e., the locking between color and flavor described by 
the group SU (2)c+l+r and the invariance under the ro- 
tations in the plane perpendicular to the external field. 
We then obtain the corresponding free energy and the 
minimum equations for the three order parameters in 
Sec. III. In Sec IV, we present our main numerical so- 
lutions for all the gaps parameters. A discussion of the 
results and our main conclusions are given in Sec. IV. 
In Appendix A, the details of the Fierz transformation 
under both Lorentz and rotational symmetry breaking is 
given. In Appendix B, the spin-1 nature of the magnetic- 
moment condensate is discussed. The Lagrangian density 
for lowest Landau levels is obtained in Appendix C using 
a chiral-spin representation. 



II. GAP STRUCTURE IN MCFL WITH 
MAGNETIC MOMENT 

If we neglect the masses of all three light quarks 
at ultra-high densities, the general CFL pairing ansatz 



Q spontaneously breaks the full QCD symmetry group 
SU(3)c x SU(3) L x SU(3) R . x U(l) B into the diago- 
nal SU(3)c+l+r subgroup. The original electromag- 
netism U(1)q is an implicit subgroup of flavor symmetry. 
Within the residual gauge group SU(3)c+l+r, there ex- 
ists a U(1)q transformation, ip — > e~ l< ^ e tp (ip is the quark 
field) that leaves the pairing invariant. The generator of 
this remanent symmetry is 

Q = Q x 1 + 1 x -Lt 8 . (1) 
v3 

with Q the conventional electromagnetic charge of 
quarks and Tg the 8th Cell-Mann matrix. Upon 
the choice of the representation of the matrices Q = 
diag(-l/3,-l/3,2/3) for (s,d,u) flavors and T 8 = 
diag(-l/ v / 3,-l/V^,2/v / 3) for (b,g,r) colors, wc find 
that the Q charges of the different quarks are [j| 
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The eigenfield corresponding to the symmetry gener- 
ated by |T]) is given by a linear combination of the photon 
field A M and the 8th gluon field 0, Q 

=cos6A IM - sin OG^, (2) 

The long-range field A^ behaves in the CS medium as the 
new effective electromagnetic field. The corresponding 
orthogonal linear combination 

&l = sin fl^ + oosflG*. (3) 

will be massive. The mixing angle is defined by cos 9 = 
g/\/e 2 /3 + g 2 , which is sufficiently small (6 ~ 1/40) at 
moderate density (g 2 /(47r) ~ 1). Effectively, the "ro- 
tated" photon ((5J) mostly consists of the usual photon 
with a small admixture of the 8th gluon. 

Clearly, the presence of an external magnetic field 
explicitly breaks the flavor symmetry from SU(3)l x 
SU(3) r to SU(2) L x SU(2) R , since the different electric 
charge of the s and d quarks from the u quark's, is man- 
ifested in the coupling with the external field. Assuming 
that the most energetically favored condensate is the one 
that preserves the highest degree of symmetry, it was pro- 
posed in Q a condensate ansatz that in the presence of 
an external magnetic field spontaneously broke the sym- 
metry to the color-flavor-locked subgroup SU(2)c+l+r- 
This implies that the CFL gap parameters Ai, A2 and 
A3, which describe d — s, u — s and u — d diquark pairs 
respectively, and which satisfy Ai = A2 = A3 in the 
CFL phase, change to A — Ai and A^ — A 2 — A 3 in 
the MCFL case. This is the MCFL ansatz Notice 
that even though all diquark pairs are neutral with re- 
spect to the rotated electromagnetic charge, they can be 
formed either by a pair of neutral or by a pair of oppo- 
site rotated charged quarks. Explicitly, the gap A gets 
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only contribution from pairs of neutral quarks, while Ab 
has contributions from both charged and neutral quarks 
and should directly feel the background held through the 
minimal coupling of the quarks in the pair with B. Al- 
though the pairs contributing to A are all neutral, this 
gap is indirectly affected by the external held through 
the coupled gap equations. The symmetric gaps are also 
separated in two in the MCFL phase, so there is one As 
and one Asb that were indeed considered in the deriva- 
tions of Ref. In the present work, however, we shall 
ignore them, as they are always smaller than their corre- 
sponding antisymmetric counterparts. 

As previously mentioned, we expect a nonzero mag- 
netic moment condensate to be present along with the 
MCFL energy gaps. That condensate can only exist 
if a proper channel opens up due to the external field. 
An external magnetic field (assumed here to be along 
the z-direction) introduces a normalized tensor — 
F^ v /\B\, with /i, v = 1,2. Because of this extra tensor 
in the theory, the metric tensor can be separated into 
transverse 

5 f = w; (4) 

and longitudinal 

9f=9^-9T (5) 

components. Next, the finite density introduces yet an- 
other normalized vector in the theory, the four-velocity 
it M , which in the rest frame reduces to = (1,0,0,0). 
With these structures, the four-fermion interaction term 
in the system Lagrangian density 

C mt = -G(#^)(#W), (6) 



with quark-gluon vertex r° = 7 M A a (where A a 's are the 
generators of SU(3)c group) , has Dirac contraction given 
by 

7 M 7 M = [au^Uv + bg^ + c{g\ v - u Al u„)] 7 'V. (7) 

Hence, one can write the four-fermion interaction ([6]), in 
the rest frame, as three distinct terms 

C mt = - 5£ (^ 7o A a V)(^7oA» - g M (^\ a ip)($lx\ a i>) 
-^(V^A^X^A'V). (8) 

Using a Fierz transformation (see Appendix A for de- 
tails) one can readily verify that the breaking of the 
Lorentz and rotational symmetries gives rise to new 
particle-particle channels of interaction, and in particular 
to (cr Qb 7 5 C)(C7 5 cr a6 ) (withC = i 7 2 7 ° anda ab = \ [l a n% 
where a,b = 1,2). Notice that the structure j 5 a ab C pre- 
serves Parity. 

What color-flavor structure can have the new con- 
densate Am — {'4' T C"f 5 a\2ip) if it has to preserve the 
color-flavor symmetry of the known MCFL ground state ? 
First, because we want to guarantee the strongest at- 
tractive channel, we choose it to be antisymmetric in 
color. Second, to ensure the total antisymmetry required 
by Pauli principle, and given the symmetric nature of 
Cj50~ab under transposition in Dirac indexes, Am should 
be symmetric in flavor. Finally, it is natural to expect 
that this condensate will receive contributions from the 
same sector of Cooper pairs that contributes to the gap 
Ab, that is, the color-flavor sector that contains pairs 
of charged quarks that minimally couple to the field in 
the Lagrangian. On the basis of all these considerations, 
we propose that the gap structure in the presence of a 
magnetic field takes the form 



I 

§«/ = Ae^%- 3 + A B (e Q/31 e u i + e a ^e ij2 ) + A M [e a01 (5 a S j3 + S l3 S j2 ) + e a ^(S a 6 j3 + fo^i)], (9) 



where a, ft and i, j denote color and flavor indices respec- 
tively. Notice that all the coefficients in ([9]) are actually 
matrices in Dirac space, for which we defined A = AC 75, 
Ab = AsC 7 5 and Am = AmC^svu- 

One can see from © that the color-flavor structure of 
the gap Am is different from that of the conventional 
gaps A and Ab, which are antisymmetric in both color 
and flavor. It is easy to check that the ansatz (|9]) is 
invariant under the same color-flavor symmetry than the 
MCFL ansatz [5[. That is, the SU(2)c+l+r symmetry, 
which requires the invariance of $y under simultaneous 
flavor (1 4-> 2) and color (lf>2) exchanges. 

The symmetric Dirac structure of Am corresponds to 
a spin-1 pairing of two quarks (ip T C'y 5 ai2tp) ■ As shown in 
Appendix B, the specific spin-1 condensate we are consid- 
ering has zero-spin projection (Ms — 0) along the field di- 



rection. Thus, it corresponds to a symmetric wave func- 
tion associated to pairs formed by quarks with opposite 
charges and opposite spins in the z-direction, and con- 
sequently with net magnetic moment projection in the 
z-direction. An order parameter with this same Dirac 
structure was considered in a 2SC model with no mag- 
netic field in Ref. [l9|, [2(| • In contrast to our case, where 
the rotational 0(3) symmetry has been already broken 
by the magnetic field, the spin-1 condensate in [19| spon- 
taneously breaks 0(3). Moreover, since only quarks of 
a single color participate in the condensate of it is 
a color-symmetric condensate which in addition breaks 
the remnant U(l) electromagnetic symmetry because has 
nonzero net rotated charge. In our case, as all the pairs 
remain neutral, the U(l) is not broken and there is no 
Meissner effect for the rotated electromagnetism. 
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We have already argued in the Introduction why a new 
order parameter Am is unavoidable on the basis of sym- 
metry arguments, because the 0(3) rotational symmetry 
is explicitly broken by the magnetic field to the 0(2) 
group of rotations about the axis parallel to the field. 
The magnetic field dependence of this condensate also 
separates it from other spin-1 condensates found in the 
literature on color superconductivity [l9|-[2l]]. In those 
cases, the spin-1 condensate is always several orders of 
magnitude smaller than the spin-0 condensate. In our 
case, as we will show below, at sufficiently high magnetic 
fields the spin-1 condensate Am is comparable in magni- 
tude to the gap Ab- 

Finally, because the Am condensate does not break 
any other symmetry that was not already broken by the 
magnetic field and by the condensates A and Ab, one 
has that the number of Goldstone fields of the MCFL 
phase [fj Q are not changed by the existence of this new 
order parameter. 



III. FREE ENERGY AND GAP EQUATIONS 

In order to obtain the free energy and to derive the gap 
equations for the three unknown order parameters A, Ab 
and Am, we use a massless three-flavor model with the 
local NJL-type interaction ©. In this model, the inter- 
action between quarks has been simplified to a pointlike 
four-fermion one, while still maintaining the symmetry of 
QCD and the same quantum numbers of the one-gluon 
exchange. The free-energy density J 7 of this system is a 
functional of the gap functions and can be obtained from 
the partition function, 



; j^ e -l3V(F+(B 2 /8tt)) 



T>i\)T>i\} exp 



d 4 x(C 



B 
8^ 



2 1 



(10) 



where N is the normalization constant, V is the 3-volume 
of the system, and (3 is the inverse absolute temperature. 
The Lagrangian density reads 



C = ^>{i$ + eQA + /i 70 )V> - (11) 

where Ap is the in-medium electromagnetic field corre- 
sponding to the symmetry U(1)q with coupling constant 
e = e cos (9, and ip is the quark spinor with implicit color 
and flavor indices. The shortcoming of this model is that 
the results will not converge in the ultraviolet region be- 
cause of the lack of asymptotic freedom, which is a sig- 
nificant asset of QCD. Hence, we have to add this feature 
by hand through introducing a UV momentum cutoff A. 
Moreover, we have another free parameter in this model, 
which is the four-fermion coupling constant G. These 
two parameters (A and G) will be fixed to reproduce the 
CFL gap Ao at a reference chemical potential for B = 0. 
In the following, we take the chemical potential and the 



momentum cutoff to be 50QMeV and lGeV, respectively. 
Accordingly, the coupling constant will be set to yield the 
CFL gap A = 25MeV. 

It is a little more involved than usual to obtain the 
particle-particle channels for the Lagrangian density (fTTj) 
via Fierz transformations. As it is known, the stan- 
dard Lorentz-covariant basis of 4 x 4 matrices is cre- 
ated from the combinations of scalar, vector, tensor, 
axial-vector and pseudo-scalar structures formed from 
the Dirac 7-matrices, i.e. {1, 7 M , <7 M ", 7^75, 175}. Now 
in a CS, Lorentz symmetry is broken down by the pres- 
ence of the Fermi sea to a mere rotational symmetry. 
Furthermore, the rotational symmetry is partially bro- 
ken by the in-medium uniform magnetic field to ro- 
tations in the plane perpendicular to the field direc- 
tion, which we assume points in the f-direction. There- 
fore, we need to work with the basis matrices given 
by { 1 , 70 , 7 Q , 7 3 , a a0 , er 30 , a ab , a 3a , 7o75 , 7°75 , 7 3 75 , «7s } 
with a, b = 1,2. This procedure will yield to (see de- 
tails in appendix A) 

£ int = G'(4>P v ^)(^P n ^) + G"^M p ^)(^M p ^), (12) 



with 



and 



(P, 



a/3 



iCj 5 e 



a/3r] 



77 = 1,2,3, 



(M r )ff = C 7 5 a ab e a(Sr (S ts S Jt + 5 it S JS ) 



(13) 



(14) 



where in (|14j) the indices r, s,t take values from 1,2,3, 
keeping a cyclic order. Notice that due to the symmetry 
breaking in the presence of an external magnetic field, 
the coefficients for both channel will be slightly different. 
However, our numerical calculations show that the results 
for all the gaps insignificantly depend on the difference 
between these two coefficients. Thus, we will consider 
G' = G" = Go in the following procedures. Next, for 
each channel we introduce a complex scalar field <f> v and 
X P respectively, and write the interaction Lagrangian as 



</>n<l>v + X*pX P 
Go 

\^P^ T )cp v + ]-^ T P v i,) 



2 

+ \(^Mp¥)x P + \x*p{VM P ^)- (15) 

Then, using a Hubbard-Stratonovich transformation, 
introducing for the complex scalar fields cj)^ and Xp 
the corresponding expectation values A, ; and (Am) p , 
and taking the mean-field approximation (i.e. neglect- 
ing the field fluctuation terms with order higher than 
two) , we obtain that the functional integral (jTU)) becomes 
quadratic in the quark fields. For the convenience of cal- 
culations, we make a rotation so (block) diagonalizing 
the gap matrix 
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which is in the basis (sb, d g , u r , s g , db, Ub, u g , s r , d r ). In- 
troducing Nambu-Gorkov spinors 



* - - ( ^(Q) \ 



(17) 



where V'(q) = ^(Q)V'j with charge projectors for the field 
representation -0* =(sb,d g ,u r ,s g ,db,Ub,u g ,s r ,d r ) given 
by 

( o) =diag(l,l,l,l,l,0,0,0,0), (18a) 
n (+) = diag(0, 0, 0, 0, 0, 1, 1, 0, 0), (18b) 
= diag(0, 0, 0, 0, 0, 0, 0, 1, 1). (18c) 

and satisfying 



(19) 



n (Q) + fi (+) + 0(_ 3 = i (20) 

the full Lagrangian density can be written as 

U , , A 2 + 2A| + 2A M 2 



Go 



(21) 

Moreover, based on the ansatz in © we already set 
Ai = A, A 2 = A 3 = A B and (A M )i = (A M ) 2 = A M , 
(Am)3 = 0. The inverse full propagator reads 

c-i m - ( [ G (Q)ol 1 $ (Q) | ( <2 9 \ 
V (Q) 1 (Q)o J / 

where Q = {0, ±} is standing for neutral, positively or 
negatively Q charged quarks. 

In (f2"2"|) we have the inverse propagators 



[Gf^r^it^^^ix-y), (23) 



with 



and off-diagonal elements 



(24) 



(25) 



with the gap matrix obtained from © after a con- 

venient rotation in color and flavor space, and given by 



(±) 



±A M - A s 

±A M - A B 



for the charged sector, and by 

,,n+ l —A 



(0) 



-A 



(26) 



(27) 



and 



(2)+ _ 

(o) - 




for the neutral sector. In the following, we will deal with 
each different block separately. The elements of the ma- 
trices (|26| - (|28|) are in color-flavor space. 

The transformation of the full quark propagator 
formed by the ©/-charged quarks to momentum space can 
be obtained by using a method originally developed for 
charged fermions in [22| and later extended to charged 
vector fields in (23|. In this approach, the diagonal- 
ization in momentum space of the Green functions of 
charged fermions in the presence of an external electro- 
magnetic field is carried out with the help of the eigen- 
function matrices TSQ^'fe), which are the wave functions 
of the asymptotic states of positive (+), and negative (— ) 
charged fermions in a uniform electromagnetic field and 
play the role in a magnetized medium of the usual plane- 
wave Fourier functions at zero field. The transformation 
functions (x) for charged quark fields are calculated 
as the solutions of the eigenvalue equation 

(n (±) . 7 )E^ ± )(x)=E^ ± )( a; )( 7 .p (±) ), (29) 

with given by 



P(±) = (po,0,±y/2eBl,p3), 



(30) 



where Z = 0, 1, 2, denotes the Landau level numbers of 
the particle with rotated charge q in the presence of the 
magnetic field 



7 ,~s G 1 

I =n- sgn(q)- + -, 



(31) 
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with n = 0, 1, 2, and a = ±1 are the spin projections. 
The EL'^ (x) functions are then given by 



E p «(or) = 4 (+) A(+) + 4- 1 <+)A(-), 



(32a) 



E#->(aO = J B p - 1 (-)A(+) + J B p (-)A(-), (32b) 

with eigenfunctions 

^(±)( x ) =7V i e-^ poa:0+J ' 2:c2+wa:3 )A(P(±)), (33) 

where -D;(yO(±)) are the parabolic cylinder functions with 
argument p defined by 



p (±) = V2lB(xi ± pa/eB), (34) 
The normalization constant Afi in (|33[) is 

A/J = (4 7 reB) 1 / 4 /Vl!, (35) 
and the spin projectors A(±) in (|32|) are defined as 

f ± 



A(±) 



(36) 



One can show that the E^- 1 functions satisfy the orthog- 
onality condition [24| 



d 4 xE l p {±) (x)E^ ± ) (x) = (27r) 4 S< 4 > (p - p')E(l), (37) 
with E' (±) = 7o (E^)t 7o , 

5^(p-p') = 6 H '5(p a -p> )5(p 2 ~p> 2 )S( P3 -p> 3 ), (38) 



and 



S(Z) = A(sgn(£B))5 10 + 1(1 - 5 a '). (39) 



Using the E^ >(x) functions, the charged fields ?A(±) 
can be transformed according to 



^(±)(«)=^^r^ c±) W(±)(p). 



"4- (27T) 4 



^ (± )(p)E p (±) ( a; ) 



(40a) 



(40b) 



where we have defined^ ^ = E~o / 

The relations (J40J) and <(29j) , as well as the orthogo- 
nality condition (J3TJ) , enable us to transform the space 
dependent part of the Lagrangian (|2ip into momentum 
space 

(41) 



where 

(42) 

The bare inverse propagators in momentum space are 
given by 

m^P'^ 1 = (27r) 4 ?W(p-y)S(0[GjJ )0 (p ((5) )]- 1 

(43) 

with [C|q)q(-P(q))] _1 formally given for neutral quarks by 



(44) 



where pm) = (poj Pi ,P2,P3) is the usual 4-momentum of 
a free particle, and for the charged quarks by 



[Gf +)0 (P(+))} 1 =p {+) ±^ 7o , 



[G(L )0 (P(-))] 1 =^ ( _ ) ±/x7 . 



(45a) 



(45b) 



with the double sign in p(±) denoting either positively 
or negatively charged particle. 

From ([43]) , we see that because of the presence of the 
operator 5(Z), the lowest Landau level (LLL) is auto- 
matically separated from the rest of the levels. In the 
MCFL phase, this separation gives rise to a spin degen- 
eracy factor for the I Landau levels, gi = 2 — Sio, in the 
free energy, reflecting the fact that higher LL's (/ > 0) 
are double degenerated. However, in the present case 
the situation is different. The magnetic field interaction 
with the pair magnetic moment breaks the spin degen- 
eracy of the nonzero LLs, a fact that will be reflected in 
the dispersion of the charged quasiparticlc excitations in 
the background of the condensate and also in the differ- 
ent way that Am enters in the off-diagonal part of the 
inverse propagator for zero and nonzero LLs. We will 
discuss more about this point in the Appendix C. 

Integrating over the Nambu-Gorkov spinors, we obtain 
the partition function as 



Q 



Det- 



VP 



1/2 



ota A 2 + 2A| + 2A m 2 , 
exp(-/3F -2 ), 

1^0 



(46) 

and hence the system free energy ft — —jjlnZ is found 
to be 



2A| +2 A M 

Gn 



(47) 



where oj k = (2k + l)n/P, k = ± 1, ±2, are the Mat- 
subara frequencies. The Matsubara sum in (j47|) can be 
evaluated using the identity 



^E ln (47^) = N + |ln(l + e-^). 



P 



IIP 2 



P 



(48) 
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In the zero-temperature limit, only the first term from 
the RHS of P5|) survives, leading to the result 



Q = - 



d 3 p 

A 2 



1 2 

5E1 



i=i 



E 

i=0 



dp 2 dp 3 
(2tt) 3 



2A| 



2A 



Go 



(49) 



where A is the energy cutoff of the NJL effective the- 
ory and rig = 7[A 2 /2eB], with /[...] denoting the integer 
part of the argument. In (I49|) . Sj and e c are the energy 
modes of the neutral and charged quasiparticles respec- 
tively, i.e., the values of the energy on which the pole of 
the determinants laid 



detS (0) [iej,p] 



detS,A [ie c 



0. 



(50) 



Notice that we use a color-flavor basis in which the gap 
matrix (j!6p is conveniently block diagonal. Therefore, 
the full inverse propagator is also block diagonal, and the 
determinant can be broken into four manageable pieces 
in (|46p. Using the identity 



det 



A B 
C D 



= det(A - B J D- 1 C)detD, (51) 



we obtain the dispersion relations of the different modes 
by finding the zeros of the determinants in (|50[) . and 
counting the corresponding degeneracy d. In this way, 
we find 



and 



ei =±V(P±A*) 2 + A 2 (d = 6), 



(52) 



e 2 =± 



2(A 



;A 2 + (^+^ 



± 1 \/ A 2 A 2 + 16[2A M V ± + m 2 (p! + *§)] ± S^A 2 ^ + A2p 2 } 



1/2 



(d = 2), 



(53) 



for neutral quarks. In (|53[) . we defined A 2 = A 2 + 8A 



,1/ 



and A 2 = A 2 +8A 2 3 . The dispersion relations for charged 
quarks in higher LL's (I > 0) are 



±\ 2eBl + A 



M 



A 2 B 



A' 2 +P3 ± 2J2£Bl(A M z + M 2 ) + (m ± A M A B ) 2 (d = 4), 



(54) 



r 



with p\ = Pi+P2- Taking into account the degeneracy d 
of each mode and the products of all sign combinations, 
we have that the total number of modes in (|52")) - ([5"4l is 
72, in agreement with the total number of degrees of 
freedom of the fields defined in color x flavor x Dirac x 
Nambu — Gorkov space. The double sign in front of Am 
in (f5"4"|) reflects the breaking of the spin degeneracy for 
the higher LL modes of the charged quasi-particles due 
to the presence of the magnetic-moment condensate Am- 
Clearly, the dispersion relation for the LLL cannot be 
found as the limit (1=0) of ([54^1. since the LLL does not 
have spin degeneracy. As in the case of massless QED 
with a dynamically generated anomalous magnetic mo- 
ment [l3l ] , the dispersion relation for the charged fermions 
in the LLL has to be found independently in a reduced 
space with only one spin projection. This is done in Ap- 
pendix C where we obtain the LLL dispersion mode 



e c = ±^/{p 3 - M ) 2 + (A B - A M ) 2 



(55) 



It can be readily checked that the dispersion relations 
(l52l - (|54l) recover the form of those previously found in 



the MCFL phase if one puts by hand A m = 0. The 
explicit breaking of the rotational symmetry by the mag- 
netic field is manifested in the separation of the parallel 
and perpendicular components of momenta in the modes 
(El and (El. 



The free energy can be explicitly found by substitut- 
ing the dispersion relations (|52l) - (j5"5)) in (|4"9"|) . A stable 
phase must minimize the free energy with respect to the 
variation of the three gap parameters, A, A# and Am- 
This gives rise to the gap equations 



an 



OA dA B OA 



M 



(56) 



These gap equations are quite complicated, even in the 
strong-magnetic-field limit where only the contribution 
from the LLL is important, and can only be solved nu- 
merically. 
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FIG. 1. Gap parameters as a function of eB/fi 2 for fi — 
500MeV. The corresponding magnetic field for eB / fi 2 — 1 is 
4 x 10 19 G. For small magnetic field below those shown in the 
figure, Am/Aq behaves as 0.32eB//i 2 . 



IV. NUMERICAL SOLUTIONS AND 
DISCUSSION 

As previously analyzed, the integral and sum in (j49|, 
and consequently in the gap equations (j5"fl)) . must be de- 
fined up to a cutoff A. However, as pointed out in in 
order to avoid unphysical discontinuities in many thermo- 
dynamical quantities, it is useful to introduce a smooth 
cutoff function h\ that should approach 1 at low ener- 
gies and at large ones. In this paper, we will follow the 
same choice as in 10] 

h A = exp(-£ 2 /A 2 ), (57) 

where for neutral quarks £ = p and for charged quarks 

£ = \J p 2 + 2eBl. This Gaussian-like cutoff implements 
the converge in the momentum integral and sum over 
Landau levels and will not leads to unphysical disconti- 
nuities. 

In Fig[TJ we plot the gaps as functions of a dimension- 
less parameter e~B/[i 2 . For small magnetic field, A and 
Ab are close to each other and approach the CFL gap 
Ao = 25MeV. As the magnetic field increases, A and 
Ab display oscillatory behaviors with respect to eB/pi 2 
as long as eB < fi 2 . These oscillations are known as the 
de Hass-van Alphen phenomenon and appear in different 
charged fermion swtems under magnetic fields (see for 
instance [l(| El, HH). As originally explained by Lan- 
dau [26j], these oscillations reveal the quantum nature of 
the interaction of the charged particles with the magnetic 



field (what is now called the Landau quantization phe- 
nomenon), and are produced by the change in the density 
of states when passing from one Landau level to another. 
The oscillations cease when the first Landau level exceeds 
the Fermi surface. For ultra-strong fields, when only the 
LLL contributes to the gap equation, Ab is much larger 
than A, which was first found by analytical calculation 
in Q. The reason for this phenomenon is that the field 
increases the density of states of the charged quarks and 
thus the pairing of charged particles will be reinforced by 
the penetrating magnetic field. 

It is apparent from the graphical representation of Am 
in FigQ] that its value remains relatively small up to 
magnetic- field values of the order of fi 2 . In the field region 
between 10 18 — 10 19 G, the magnitude of Am grows from 
a few tenths of Mev to tens of Mev. It becomes compara- 
ble to the MCFL gap A b when the field is strong enough 
to put all the quarks in the LLL, shown in the final seg- 
ment of the plots in the figure. The de Hass-van Alphen 
oscillations are much smaller for Am than for the regular 
gaps. These features indicate that the main contribution 
to this gap should come from pairs whose charged quarks 
are at the LLL. From a physical point of view, this can 
be understood taking into account that the contribution 
from higher LL's (I > 0) should be negligible because the 
magnetic moment of a pair where the positive quark has 
spin up and the negative quark has spin down cancels 
out with that of a pair where the positive quark has spin 
down and the negative quark has spin up. The only con- 
tribution from higher LLs can come when the number of 
particles is odd, so there are energy states occupied by 
a single particle, but that is a very small part. The can- 
celation does not occur, however, between the pairs of 
quarks in the LLL because they can only be formed by 
positive quarks with spin up and negative quarks with 
spin down. At low fields, the number of quarks in the 
LLL is scarce, while for fields of order eB > /i 2 , all the 
particles are constrained to the LLL, hence the variation 
of Am from lower values at weak field, to higher values 
at sufficiently strong fields. 

Another important consequence of the new gap is the 
increment in the magnitude of As for any given value of 
the magnetic field in the strong field region, as compared 
to its own value found at the same field but ignoring the 
existence of Am- This effect, combined with the increase 
of Am at strong fields, will make the CFL phase in the 
presence of a strong magnetic field more stable than the 
regular CFL, a fact that could favor the realization of an 
MCFL core in magnetars. 

Our calculations show that at zero magnetic field Am 
is zero, so the system does not behave as a ferromagnet. 
This result is consistent with the rotational invariancc 
of the CFL phase because although the Cooper pairs 
formed by charged quarks have nonzero magnetic mo- 
ment even in the absence of a magnetic field, if these 
moments were spontaneously aligned in some direction, 
they would break the rotational symmetry of the CFL 
phase. Accordingly, the expectation value of the mag- 
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FIG. 2. The ratio 47reA///i 2 as a function of eB/^i 2 with M 
the magnetization for /j, = 500MeV and A m 7^ in the solid 
line; and with Am = in the dashed line. 

netic moment at zero field (in the CFL phase) must van- 
ish: 

(0\A m \0)cfl = (OIR^RAmR^R^cfl 

= (Q\RA m R- 1 \0)cfl = 0. (58) 

However, the existence of expectation value of Am in the 
presence of a magnetic field is unavoidable. We found 
that there is no solution of the gap equations with B ^ 
0, A ^ 0, A B ^ 0, and A M = 0. 

The difference found between the A's and Am gaps is 
also in agreement with the fact that the extra order pa- 
rameter Am is a consequence of the magnetic field (i.e. 
Am = at B = 0). In contrast, the relevant scale for 
the generation of the A's gaps is the energy at the Fermi 
surface, i.e. the chemical potential. It is logical that once 
the magnitude of the magnetic field is comparable to the 
chemical potential, the induced magnetic moment con- 



If the value of the magnetic field is such that the Fermi 
sphere can accommodate n LLs, the dispersion of the 
last level occupied within the Fermi surface can be either 
fcrmionic or bosonic, depending on whether the level is a 
little below or exactly at the Fermi surface. For example, 
in the case that the Fermi sphere can exactly accommo- 
date two LL's (i.e. for field values eLB//i 2 ~ 0.5), the level 
/ = 1 lies at the Fermi surface. Under this condition, the 



densate becomes as large as the gap. Another important 
consequence of the generation of the magnetic moment 
condensate is that its presence strengthens the gap Ab 
in the sufficiently strong-magnetic-field region, as can be 
checked by comparing our results in Fig. 1 with those of 
Ref. [H. 

Finally, the magnetization M = —dT/dB, which 
describes the average magnetic moment per unit volume 
of the superconducting medium at the stationary point, 
is depicted in Fig. [2] as a function of the magnetic field. 
Comparing the two cases, Am = (dashed line in Fig. 
[2]) and Am ^ (solid line in Fig. [5]), it is evident 
that at strong fields the magnetization is reinforced by 
the magnetic moment of the Cooper pairs. It increases 
about 10%-20% for fields in the range fi 2 < eB < 2/i 2 , 
as shown in the shift between the two curves in the 
figure. The magnetization also exhibits, like the gaps, 
an oscillatory behavior. This larger magnetization at 
strong fields will be reflected in the equations of state of 
the magnetized system through the transverse pressure 
1. 



V. DISPERSION RELATIONS 

Let us discuss the effect of the magnetic field on the 
spectrum of the charged quasiparticles. In Fig. O we 
plot the dispersion relation of a quasiparticle in the LLL 
(positive energy in Eq. ([55]) ). At the field value con- 
sidered in the figure, eB/[i 2 ~ 2, the LLL is the only 
level that fits in the Fermi sphere. As can be seen there, 
the spectrum of the LLL corresponds to that of the BCS 
phase with a minimum at p^ — fi. 

For higher LL's we need to use the dispersions given 
in Eq. (|54|) . There, the overall double sign denotes par- 
ticle/hole around the Fermi surface, the double sign in 
front of the inner square root is for particle/antiparticle, 
and the double sign in front of Am corresponds to the 
two possible projection of the spin of the quasiparticle in 
the magnetic field. In the case of a quasiparticle, with 
magnetic moment opposite to the field we have 



(59) 

I 

quasiparticles at 1=1 has minimum energy at p-j = 0, a 
bosonic type of dispersion as shown in Fig. 2J If we de- 
crease the field just a little, 2eB < [i 2 , still there would be 
only two LLs inside the sphere, but the minimum of the 
dispersion for I — 1 would be at a nonzero p^, and thus 
/ = 1 would behave as a fermionic mode. If we continue 
decreasing the field to eBf/j 2 ~ 0.25 we find the same sit- 
uation, that is, that the upper LL available in the Fermi 



I 

2eBl + A 2 M + A| + fi 2 + p\ - 2\fle~Bl{A 2 M + y 2 ) + (/ip 3 + A m Ab) 2 I > 1- 
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FIG. 3. Dispersion relation for charged quarks at eB/fi 2 = 2. 
At that field value the quasiparticles are in the LLL and show 
a BCS spectrum. 



sphere (the one with I — 2 in this case) will have a bosonic 
spectrum, while the other two have BCS spectra, as can 
be seen in Fig. [5] It is easy to understand that the same 
behavior will be found at any field value. The reason is 
that once the Fermi sphere is filled with more than one 
LL, if the upper LL is already at the Fermi surface, then 
it is only needed an energy equal to the gap to excite a 
quasiparticle, and consequently the minimum of the dis- 
persion relation takes place at p% — 0. The exception 
occurs when only the LLL fits in the Fermi sphere. In 

this case, there is no transverse momentum pj_ = VzeEBZ 
to equate the Fermi sphere radius \i since 1 = 0. Hence, 
to excite the LLL quasiparticle from the Fermi surface a 
longitudinal momentum equal to the Fermi sphere radius 
P3 = fi is needed, plus the energy gap. 

As known, the analytic behavior of the quasiparticle 
spectrum can determine the nature of the diquark con- 
densate [13, . It turns out that a magnetic field mod- 
ifies the quasiparticle dispersions and consequently can 
change the nature of the diquark pairs. In [29j the effect 
of a magnetic field in the dispersions of the quasiparti- 
cles in superconductivity was found to be equivalent to 
producing an effective mass Mi = 2eBl for the quasipar- 
ticles. The relation between that effective mass and the 
chemical potential determines the fermionic (p > Mi) or 
bosonic (/i < Mi) character of the mode. As discussed 
in |2!jj . a magnetic field can tune the crossover between 
a BEC and a BCS regime in a color superconductor by 
varying the relative numbers of LLs for which the effec- 
tive chemical potential fj,i = fj, — Mi is either positive or 
negative. 

An important outcome of our calculations is the en- 
hancement of the net energy gap of the LLL quasipar- 
ticles at strong magnetic fields. As can be seen from 
([55]). the energy gap of the LLL quasiparticles is given 
by Ab — Am- Fig Q] shows that the condensate Am is 
negative for all field values, has a small magnitude from 
low to intermediate fields, and then becomes comparable 
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FIG. 4. Dispersion relations for charged quarks at eB/fi 2 = 
0.5. Different lines correspond to the available LL's I at that 
field value, with black-solid line for I — 0, and red-dashed 
line for I = 1. The upper available LL (I — 1) has a BEC 
spectrum. 



in magnitude to in the region of strong fields. On the 
other hand, as mentioned above, the gap Ab, which was 
found to be enhanced by the magnetic field in the strong 
field region in the case where Am was taken zero by hand, 
becomes even larger in the presence of Am- Therefore, 
the net energy gap at strong fields is much larger than 
at zero field. This enhancement effect means that the 
condensation energy of the pairs of charged quarks is en- 
hanced and this in turn can have implications for the 
chromomagnetic instability associated to the existence 
of gapless modes at moderate densities when the mass of 
the s-quark, M s cannot be neglected and the color and 
neutrality conditions are enforced [13, 30]. The criterion 
of the phase transition from the CFL to the gapless-CFL 
phase was found to be [30j 

M 2 

— - « 2A CFL . (60) 
A* 

If the gap on the RHS of is enhanced through any 
effect (in our case it would be due to the presence of 
a strong magnetic field), the CFL (MCFL if a field is 
present) phase can bear larger values of M s without be- 
coming unstable, thus delaying the transition to a gapless 
phase to even lower baryon densities. Even though the 
model we considered in our work ignores the effects of 
the strange quark mass, the arguments used to justify 
the existence of an extra condensate Am in the presence 
of a magnetic field will not change if M s is included. 

One can check that the spectrum (15"3l) . associated to 
the neutral quasiparticles that couple with the charged 
ones through the gap equation, has a larger energy gap 
due to a nonzero Am- Hence, at strong fields, Am also 
produces an enhancement in the energy gap of some of 
the neutral excitations. We expect that this effect will de- 
lay the appearance of gapless excitations for those modes 
too when the strange mass is included in the analysis. 
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FIG. 5. Dispersion relations for charged quarks at eB//j, 2 = 
0.25. Different lines correspond to different available Lan- 
dau levels I at that field value. Here, the different spectrum 
correspond to: black-solid line (1 = 0), blue-dash-dotted line 
(I = I), and red-dashed (1 = 2), respectively. Notice that only 
the upper available LL has a BEC spectrum. 



VI. CONCLUDING REMARKS AND OUTLOOK 

In this paper, we extend the study of the MCFL phase 
of color superconductivity to include the effect of the 
magnetic moment of the Cooper pairs formed by rotated 
charged quarks. We found that the external magnetic 
field aligns the magnetic moments of those pairs produc- 
ing a net expectation value of the magnetic-moment of 
the pairs, Am/\B\, in the field direction. The new order 
parameter A M leads to various effects. First, the gap An 
is larger in the strong field region than what it was in [5| , 
where Am was not considered. In addition, the net en- 
ergy gap of the charged quasiparticles in the strong field 
region is more than twice the one in |5[ . Finally, there is 
no solution that minimizes the MCFL free energy with 
A ^ 0, A B ^ and A M = 0. Thus, in the MCFL phase, 
the magnetic moment condensate has to be considered in 
equal footing with the gaps A and A b ■ 

The justification for the increase in the number of or- 
der parameters in the MCFL phase is easy to understand 
on the basis of symmetry arguments. The reason is that 
the magnetic moment condensate Am does not break 
any additional symmetry that has not already been bro- 
ken by the gaps A and A# and the magnetic field, thus 
there is no reason why Am needs to be zero. The sit- 
uation has some resemblance with the phenomenon of 
magnetic catalysis in massless QED, where an applied 
magnetic field explicitly breaks the rotational symme- 
try and also catalyzes the spontaneous breaking of chiral 
symmetry through the generation of a chiral condensate. 
Once these two symmetries are broken there is no reason 
why any physical parameter whose presence in the La- 
grangian would break the exact same symmetries has to 
be forbidden. Consequently, two parameters, the mass 
and the magnetic moment are dynamically generated in 
the QED case [Hj]. 



Technically, the emergence of a new condensate in the 
MCFL phase is connected to the modification of the 
Fierz identities in the presence of a magnetic field. To 
see this one should first notice that the Lorentz sym- 
metry is broken by the dense medium down to merely 
a rotational symmetry and, in addition, the magnetic 
field further breaks this symmetry from 0(3) to 0(2). 
When the Fierz transformations are performed in the 
point-like four-fermion interaction of a NJL theory with 
these explicitly broken symmetries, one immediately find 
the opening of new pairing channels, one of which fa- 
vors the formation of a spin-1 condensate of the form 
Am = (V ;T CS 3 7 5 i/)), which, as previously discussed, is 
unavoidable. 

We call the reader's attention to the fact that the 
Am condensate of the present work is quite different 
from other spin-1 condensates previously considered in 
color superconductivity. Usually, spin-1 condensates not 
only are much smaller than the s-wave ones (2lj . but 
they do not mix up with s-wave gaps in the quasiparti- 
cle spectrum. The existence of quasiparticle dispersions 
with very small gaps in other phases with spin-1 conden- 
sates makes them relevant for the transport properties 
of the system, but this same characteristic makes these 
condensates very sensitive to be erased by temperature 
and, in cases where they break the rotated electromag- 
netism, also by magnetic fields. In contrast, the spin-1 
condensate of our case mix up with the s-wave gaps in 
the quasiparticle spectrum in such a way that at strong 
fields it actually increases the effective energy gap of the 
quasiparticles. Besides, the extra condensate has also 
the indirect effect of making the gap A b even larger in 
the strong field region, hence raising the critical tem- 
perature required to erase the color superconductivity. 
This feature makes this phase of potential interest for 
the planned low-temperature/high-density ion-collision 
experiments at different worldwide facilities, as pointed 
out in the Introduction. 

In our work we neglected the strange quark mass M s . 
However, this may not be an accurate assumption for 
the intermediate densities prevailing in the interior of 
neutron stars. The response of the CFL superconduc- 
tor to a nonzero strange quark mass and the constraints 
imposed by color and electrical neutralities can lead to 
gapless dispersion relations when the density decreases 
[30l | . This in turn has been found to produce chromo- 
magnetic instabilities [17j. We do not see any obvious 
reason to expect that the enhancement of the energy gap 
at strong fields found in the present work will not be 
present too if the strange quark mass is taken into con- 
sideration. If the energy gap of some of the quasiparticle 
dispersions is much larger, obviously the density needed 
to turn them into gapless modes (and hence produce the 
chromomagnetic instability) will be much smaller, a fact 
that could favor the realization of the MCFL phase in the 
core of a strongly magnetized neutron star, even at inter- 
mediate densities. Therefore, an important pending task 
will be to investigate the MCFL considering a nonzero 
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M s and imposing the neutrality conditions, while taking 
into consideration the effect of the Cooper pairs' mag- 
netic moment. 

Even if the MCFL phase does not turn out to be the 
most favored phase in neutron star cores or in the planned 
experiments, the present investigation is also important 
from a fundamental point of view, because it has uncov- 
ered a mechanism to increase the energy gap that could 
be also applicable to other color superconducting phases 
more relevant at intermediate densities. For example, it 
would be interesting to explore whether a spin-1 conden- 
sate similar to the one found here is also favored in the 
strongly coupled 2SC phase in the presence of a magnetic 
field. It has been argued that the 2SC phase at strong 
coupling can be the most energetically favored at mod- 
erate densities [3l|. All the Cooper pairs in this phase 
are formed by quarks with opposite spins and charges, 
so they all have nonzero magnetic moments which pre- 
sumably could be oriented by an external magnetic field. 
This orientation could lead to a net magnetic moment, 
likely favoring the formation of a new condensate, as in 
the MCFL case. If, as in the MCFL situation, the extra 
condensate in the strongly coupled 2SC theory makes the 
regular s-wave gap even larger, this will lead to a larger 
critical temperature, making this phase a potentially in- 
teresting candidate for the QCD phases that will be ex- 
plored by the planned low-temperature/high-density ex- 
periments. 

From Fig. [2]we saw that the magnetization is enhanced 
at strong fields. This feature could serve to determine 
whether the core of magnetars is made of color super- 
conducting matter or of hadronic matter, whose magne- 
tization is known to be negligible even at large fields [32| ■ 
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Appendix A: The Fierz Identities with both Lorentz 
Symmetry and Rotational Symmetry Broken 

In this Appendix, we present the details of the Fierz 
transformation with both Lorentz and rotational symme- 
try breaking. Fierz identities are connected to reordering 
of field operators in a contact four-particle interaction. 
Let us consider an interaction: 

(^A^X^-BM. (Al) 

Here, the indices of the spinors are suppressed. The 
same interaction can be expressed in a different way 



as (i/iiM^Jf^Af^)- How the matrices A, B relate to 
M, N is the general aim of Fierz transformations. 

The 16 Dirac bilinears in (|A1[) are usually classified 
into distinct classes according to their properties under 
Lorentz transformation as 

{r^} = {1,75,7^,757^,^}, (A2) 

The orthogonality relation for the basis {T A } are 

Tr[r A r B ]=4^, (A3) 

This relation allows us to expand any complex 4x4 ma- 
trix X in terms of the basis (|A2[) as 

X = X A T A , X A = ^Tr[XT A ], (A4) 

Combining (|A3[) and (IA4|) . extracting each element of the 
matrix, we could find a completeness relation as 

^mi^nj = ^(r j 4)nm(r (A5) 

This identity is sufficient to reproduce all possible Fierz 
identities by appropriately incorporating identity matri- 
ces. For general matrices X,Y, we have 

X tj Y kl = (X%(iy)ju = X -(XY c Y) a {Y c ), k 

= ^Tr[xr c Yr D ](r D ur% k , (A6) 

In particular, if X = T A and Y = T B , Eq. ((Ml) will lead 
to the Fierz identities [33| 

(r A Mr B ) fc; = ^TT[T A T c r B r D }(T D ur c ) Jk . (A?) 

Where, all the lower case letters k, I run over 0, 1, 2, 3. 
The only remaining work is to calculate the expan- 
sion coefficients which are straightforwardly obtained as 
gamma matrix traces. For the particle-particle channel 
(-0MCi/j T )(^ T CA''i/)) as in CS, the Fierz identity are re- 
lated to the coefficients in (|A7I) simply by the sign, which 
defines the (anti)symmetry of the representation CT A un- 
der transposition (when A = B). 

For the CFL phase, we should take into account 
that the Lorentz symmetry is broken by the presence 
of the dense medium down to mere rotation symme- 
try. One thus needs to work with the basis of rotation- 
covariant matrices, { 1 , 70 , 7° , <r a0 , cr ab , 7075 , 7^5 , i-f 5 } , 
with indices a and b running from one to three. 
Moreover, for the MCFL phase the rotation symme- 
try will be further broken by the penetrating mag- 
netic field that is taken along the z-direction. There- 
fore, the expansion basis (|A2[) has to be changed to 
{ 1 , 70 , 7 a , 7 3 ,cr a0 , cr 30 , cr ab , cr 3a , 7075 , 7 a 75 , 7 3 75 , 75 } with 
a,b= 1,2. 

Following the same procedure as that in going from 
(|A2|) to (|A7|) . we can obtain the Fierz identities for 
particle-antiparticle channel as 
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V («75)iz(«75)fcj / 



(A8) 



and that for particle-particle channel as 
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(a 30 ) 
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(7075) 
(7 a 75) 
(7 3 75) 

(*75) 
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i 3 {lal5)kl 
3. 



y(7 75)fc! 
»j(*7s)fc! 



)fc! / 



/ -- -- 
/ f f 

I J 



! ! 



'! 
4 

1 



1 _ _L 

I _! 

4 4 

I ? 



1 _1 _1 1 

'! f ! ! 
! f 1 ! 



? _! 4 n f 
I _l 1 _! 
I I I I 



f f 



i f 



4 

1 

I 



1 ] 

I J f 

"0 -I I 

1 ? ? 



4 4 

1 f 
I I 



4 4 



-I -t 

1 ? ? 



1 1 



1 f 



ittitnt-tntnit 

I l 1 _i 1 _| f 1 f 1 _f ? 

t J J J I I I I I I f _f 

444444444444 
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(a 30 C) 
i((i Qb C) 
(a 3a C) 
(7075C) 
(7°75C) 
(7 3 75C) 
V (*7 5 C) 



(C7o)fej 

(<?7a)fcj 

(C7s)fei 

(Cc a o)fcj 
(C<T3o)fej 
il(C<J a b)kj 

(Ca 3a )kj 
(C7o7 5 )fcj 
(C7a7 5 )fcj 
(C7 3 7 5 )fej 



(A9) 



In the following, we will start from the color current 
interaction to obtain (IT21) . Considering the Lagrangian 
density © 



£ = - 3b(V'7oAqV')(V'7oAq'(/') - 3m(V'7 J " A V)(V>7-J-A a ^) 
- ff 3 f (^ 7 3 Aa^)(^73A a ^), (A10) 



(7 X )(7±) = 5 {(<T 30 C)(Ca 30 ) + (*75C)(zC 75 ) + -} , 



1 



(Allb) 



(7 3 )(73) = — A {(a 30 C)(Ca 30 ) - (ij 5 C)(iC l5 ) + ...} , 

(Allc) 

where, the Dirac indices have been suppressed. The Fierz 
identities for the generators of SU(N) can be found in 

For particle-particle channels, we have 



which is in agreement with the symmetries if the MCFL 
phase, and taking into account the Fierz identities for 
Dirac matrices that can be extracted from (IA9I) as 



(7o)(7o) - t {(a 30 C)(Ca 30 ) - (i 7 5C)(iC 75 ) + -} 



(Alia) 



N 



N 2 -l 



2 \ f (As)ap (As) T /3 



(l)pr \ = 

(A a ) Q/ 3 (A a ) pr J 

(A12) 

Here, the Greek letters run from one to three and 
As, Xa are the symmetric and antisymmetric generators 
of SU(N) respectively. Taking into account everything 
together, the 4-fermion Lagrangian density (|A10[) be- 
comes 



C = -f-P [G'^a 30 CXs\A^ T )(^ T Ca 30 \ s \ A i,) + G" ' {^ l5 CX A X A 4 T ){^ \C l5 X A X A ^)} , (A13) 
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where we have defined the new coefficients G = \qe — 

9m + yii and G " = 9e + ^9 M + 9m- Incorporating 
the color-flavor ansatz we proposed in Section II, we can 
obtain the Lagrangian density we used in (|12[) . 



Appendix B: Spin-1 Magnetic-Moment Condensate 

Fermion condensates in relativistic theories have to 
preserve the Lorentz's covariance of the Lagrangian den- 
sity. For the chiral condensate, for example, the allowed 
Dirac structures entering in the particle-antiparticle bi- 
linear are the 16 elements of the so-called Dirac ring 

< i> a M a ^fi >, M a0 = M S I + M%>f + 

+ M* l5 Y l + M p i lb (Bl) 

Here, the supra-indices in the expansion coefficients, 
S, V, T, A and P, stand for the scalar, vector, tensor, ax- 
ial vector and pseudoscalar nature of the structure, re- 
spectively. An essential property of this decomposition is 
that each element transforms under a Lorentz transfor- 
mation into itself. For the superconducting Cooper pairs, 
the analogous Lorentz decomposition that preserves the 
relativistic covariance is 

< iZ&ccpjp > , A aP = [A S I + + 

+A£>"" + + A p J75 ] ( 75 C) (B2) 

The magnetic moment condensate corresponds to the 
tensorial element, a^ v , of the Dirac ring in (|B2[) . Since 
this condensate is symmetric in Dirac, it has wave func- 
tions \S, Ms), with S and Ms representing the pairs' spin 
and spin projection along the z-axis respectively, that can 
be given in term of the quarks' S* z -spin projection repre- 
sentation, \mi,m,2) as [34j 



(if; 1 A(+)C£VA(+)V) + W A(-)C£VAR</;), 



IM) = I+, 



11,-1) 



(B3) 



The pairs with |+,+) corresponds to quarks with both 
spins up (+), | — , — ) with both spins down (-), and |±, =p) 
to those with one quark with spin up and the other down, 
and viceversa. 

The spin-1 magnetic- moment condensates in general 
can be given by 



where E 



(V/CSVV), k= 1,2,3, 
\e k% ^<Jij is the spin operator. 



(B4) 



Introducing the sum of the spin projectors A(±) in 



[A(+) + A(-)]C£ fc 7 5 [A(+) + A(-)ty), (B5) 

and taking into account that A(±)A(±) = A(±), 
A(±)A(=f) = 0, we see from (|B5|) that for each spin 
index k, it is obtained 



A(+)CE i 7 5 A(+)V') + {V A(-)CE i 7 5 A(-)^), 



(B7) 

(V> T A(+)CE 3 7 5 A(-)</>> + (^ T A(-)CE 3 7 5 A(+)</>>, 

(B8) 

Because 

CE X 7 5 = A(-)-A(+), CE 2 7 5 = -^ (B9) 
we finally find that 



|i, i) = -^(v/ceV^) + ^ T C£ 2 rV> 



(BIO) 



|1,-1} = i^CES 5 ^) + ^ T CE 2 7 5 ^) 



= Wf-)V'(-)} 



(Bll) 



|1,0) = ^ T A(+)CE 3 7 5 A(-)V) 
+ C0 T A(-)CE 3 7 5 A(+)V^) 



<vi } 7 V (+) ) 



(B12) 



where ip^) represents the wave function of the spin- 
up (+) quarks and spin-down (-) quarks, and — 
tp T A(±)C. We can see that to get the condensate pro- 
jections Mg = 1,-1, we need to consider linear combi- 
nations of the magnetic-moment condensates along the x 
and y axes, as one would expect, since the spin operators 
in different directions do not commute. The magnetic- 
moment condensate Am we are considering in this paper 
corresponds to the order parameter (|B12[) . and conse- 
quently to a spin-1 condensate with zero component of 
the spin |1,0) in the field direction. 



Appendix C: The Lowest-Landau-Level Lagrangian 

As we mentioned above, the inverse propagator for the 
LLL should be different from those for the remaining 
LL's. In this Appendix, we will justify that statement 
and obtain the dispersion relations for the quasiparticles 
in the LLL. Considering now only the contribution of the 
positively charged quarks, we have 



£ = J d 4 x^ (+) (x)S- + )(x)^ (+) (x), (CI) 

where the NG spinors and the inverse propagator are 
given in (|17p and (|2"2"1) . Expanding this Lagrangian in 
NG space, we have 



£ 



d A x 



ip(+)(i$ + eA^ + M7o)i/>(+) 



(B6) 



+ l/>C(-)$~t + )ll>(+) + c.c 



(C2) 
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Using Ritus' transformation to momentum space for 
charged fields 



(C3b) 



and taking into account (J2H]) and (13"T)) . we have for the 
free-propagator part in (|C2I) 



7^+) + M7o V>(+), (C4) 



where S(Z) was defined in (|39|) . as the operator that sep- 
arates the LLL Lagrangian from the rest 



.,.! / d p dp 2 dp 3 

C F =C F + ^j (2?r)4 V(+) 



/=i 



7/iP( +) + M7o 



(C5) 



with the LLL part given by 

-■0(+) A (+)[7"P|| +M7o]V'(+)- (C6) 



£ o _ / dpodp2dp3 



(2tt)' 



The spin projectors A(±) can be expressed in terms of 
Pauli matrices as 



where 



A(±) = 



^ = £(1 + ^) 



With the chiral projection operators 



R = 



1 + 75 



L = 



1-75 



(C7) 



(C8) 



(C9) 



in the chiral representation, 75 = diag{— 1, 1}, we have 
the following commutation relations 



[A(±),i] = [A(±) fJ R]=0. 



(CIO) 



In terms of these projectors, we can represent the Dirac 
spinor as 



^ = v4 +) +^ ) +v4 +) +</>ir\ 



(Gil) 



with 



= flA(±ty, Vi ±} - LA(±)^. (C12) 



and then the Lagrangian density of the LLL (IC6[) can be 
written as 



dp a dp 2 dp, 



(27T) 



gr{^( + )]k +) (7 ll P||+M7o)[^ (+ )] 



(+) 



+ ^( + )]i +) (7 ll P||+M7o)[^( + )]i +) } 



(C13) 



Here, we have used the relations 

7 llA(±) = A(±) 7 H, 7 ± A(±)=A(t)7 ± , (G14) 

In terms of the explicit components of Dirac spinor 
wf+\ — ("01; i } 2, i/>3, 1P4), we can rewrite (|C14|) as 

f dp dp 2 dp 3 



(27T) 



■^(Po+M-Ps)^]- (C15) 



For the term involving the gaps in (|C2I) . we have 
d A p 



^(2^4^c(-)S(/)$+ +) V( + ), 



(C16) 



From where we can isolate the LLL part 

f dp dp 2 dp 3 T A ,,s^+ , 
L "= ' ( 2 7t) 4 ^c(-)A(+)^ +) V(+) 



dp dp 2 dp 3 
(2tt)* 



{fe H ]i i +) A(+)$+ +) ^ (+) ](; 



(+) 



+ [^ ( _ ) ]i +) A(+)$+ +) [V, (+) ]i +) }. (C17) 

where we used that A(+)A(+) = A(+). 

Notice that the gap matrix ^J.) can be written as 

= (-A B + A M )7sA(+) + (-A B - Am)7sA(-), 

(C18) 

Then in (ICTTl) only the first term in the RHS of ([Ul8| 
can survive since A(+)A(— ) = 0. Thus, we have 

£ o = f ^ ^ 2 ^ 3 [^(-A b +A m )V'3+^(-A b +Am)V'i] 

(C19) 



(27T) 



Here, ^ 4 are the components of spinor V'(-)- The same 
procedure can be applied to the charge-conjugate part in 
(|C2|) . Putting everything together, we have 



£° = C F +£° G +c.c 
with 



dp dp 2 dp 3 



*o(p)[So(p)]- 1 *o(i>), 
(C20) 



*o(p) 



^3 



(C21) 



and 



[5o(p)]- 1 -7o 



(p + A*) + 0"3P3 -A B + A M 

-A B + A M bo - M) - 03P3 



(C22) 

Now we could readily obtain the dispersion relations for 
LLL through det[5o(p)] _1 = as 



el = ±V(P3 - M) 2 + (Ab - A M ) 2 



(C23) 
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For the negatively charged quarks, the dispersion rela- 
tions are the same since the sign in (|3T)| will change ac- 
cordingly. 

Note that in the LLL Lagrangian (jC23[) the spinor 
(|C21|) is a 4-component field. This is different from the 



LLL Lagrangian in QED [13(, where the LLL fermions 
behaves as the two-component spinors of the (1+1)-D 
Thirring model (35j |. This difference is due to the con- 
tribution in the inverse propagator (|C22|) of the Nambu- 
Gorkov off-diagonal fields corresponding to the charge 
conjugate of the spin-up negatively charged quarks. 
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